o 
o 
o 

o 

Q 



A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

02(02.01.3), 03(03.20.8), 08(08.09.2 Capella,08.09.2 Procyon, 08.03.5, 
08.12.1, 08.01.2), 13(13.25.5) 



ASTRONOMY 

AND 
ASTROPHYSICS 



Helium-like triplet density diagnostics: 

Applications to CHANDRA-LETGS X-ray observations of Capella and 
Procyon 

J.-U. Ness 1 , R. Mewe 2 , J.H.M.M. Schmitt 1 , A.J.J. Raassen 2 3 , D. Porquet 4 , J.S. Kaastra 2 , R.L.J, van 
der Meer 2 , V. Burwitz 5 , and P. Predehl 5 

Universitat Hamburg, Gojenbergsweg 112, D-21029 Hamburg, Germany 

2 Space Research Organization Netherlands (SRON), Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands 

3 Astronomical Institute "Anton Pannekoek", Kruislaan 403, 1098 SJ Amsterdam, The Netherlands 

4 CEA/DSM/DAPNIA, Service d'Astrophysique, CEA Saclay, F-91191 Gif sur Yvette Cedex, France 

5 Max-Planck-Institut fur Extraterrestrische Physik (MPE), Postfach 1603, D-85740 Garching, Germany 

received February 1, 2008; accepted December 01, 2000 



> 

m : 

(N , 

(N 

(N 

o : 
o ■ 

Oh 1 
I 1 

o 1 

a: 



Abstract. Electron density diagnostics based on the 
triplets of Helium-like C v, N vi, and O vn are applied 
to the X-ray spectra of Capella and Procyon mea- 
sured with the Low Energy Transmission Grating Spec- 
trometer (LETGS) on board the Chandra X-ray Ob- 
servatory. New theoretical models for the calculation of 
the line ratios between the forbidden (f), intercombina- 
tion (i), and the resonance (r) lines of the helium- like 
triplets are used. The (logarithmic) electron densities (in 
cgs units) derived from the f/i ratios for Capella are 
< 9.38 cm" 3 for Ovn (2a upper limit) (f/i=4.0±0.25), 
9.86 ± 0.12 cm" 3 for Nvi (f/i=1.78±0.25), and 9.42 ± 
0.21 cm" 3 for Cv (f/i=1.48±0.34), while for Procyon we 
obtain 9.28^j 8 cm" 3 for Ovn (f/i=3.28±0.3), 9.96 ± 
0.23 cm" 3 for Nvi (f/i=1.33±0.28), and < 8.92 cm" 3 for 
Cv (f/i=0.48±0.12). These densities are quite typical of 
densities found in the solar active regions, and also pres- 
sures and temperatures in Procyon's and Capella's corona 
at a level of T ~ 10 6 K are quite similar. We find no ev- 
idence for densities as high as measured in solar flares. 
Comparison of our Capella and Procyon measurements 
with the Sun shows little difference in the physical prop- 
erties of the layers producing the C v, N VI, and O VII emis- 
sion. Assuming the X-ray emitting plasma to be confined 
in magnetic loops, we obtain typical loop length scales of 
Lcapeiia > 8 L Procyon from the loop scaling laws, imply- 
ing that the magnetic structures in Procyon and Capella 
are quite different. The total mean surface fluxes emitted 
in the helium- and hydrogen-like ions are quite similar for 
Capella and Procyon, but exceed typical solar values by 
one order of magnitude. We thus conclude that Procyon's 
and Capella's coronal filling factors are larger than corre- 
sponding solar values. 
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1. Introduction 

The hot plasma in the corona of the Sun and of other 
stars is thought to be in "coronal equilibrium". Atomic 
excitations occur through collisions with electrons. The 
excited atoms decay radiatively and the emitted radiation 
escapes without any further interaction with the emitting 
plasma. As a consequence, the emission is optically thin, 
and the total flux emitted in some spectral band or in a 
given emission line is proportional to the emission measure 
EM, defined as the integral of the square of the plasma 
density n over the emitting volume elements dV through 
EM — Jn 2 dV. Thus, observationally, the contributions 
of density and volume to a given observed value of EM 
cannot be disentangled. 

Stellar X-ray surveys carried out with the Einstein and 
ROSAT satellites have shown an enormous range of X- 
ray luminosity (Lx) for stars of given spectral type (cf., 
Vaiana et al. (1981)| , |Schmitt (1997)D - Typically, one ob- 
serves star to star variations in Lx of up to four orders 
of magnitude, with the largest X-ray luminosities found 
among the stars with the largest rotation rates. While 
one definitely finds a correlation between mean coronal 
temperature and X-ray luminosity ( {Bchmitt ct al. (1985) 



[Schmitt (1997) ), it is also clear that the single- most im 
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portant factor contributing to the large variations in Lx is 
the variation in emission measure. The conclusion there- 
fore is that active stars (can) have a couple of orders of 
magnitude higher coronal emission measure, while main- 
taining the same optical output as low-activity stars like 
our Sun. 
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The emission measure is directly linked to the struc- 
ture of stellar coronae, if we assume, going along with the 
solar analogy, that the X-ray emitting plasma of a stellar 
corona is confined in magnetic loops. The observed values 
of EM and Lx for a given star could be accounted for ei- 
ther by the existence of more loops than typically visible 
on the solar surface, by higher density loops or by longer, 
more voluminous loops. Thus the question is reduced to 
the following: if EM* 3> EMq for an active star, one 
wants to know whether this is due to n± > n Q or T4 > Vq 
or both. 

Spatially resolved solar observations allow to disen- 
tangle density and volume contributions to the overall 
emission measure. One finds the total X-ray output of 
the Sun dominated - at least under maximum condi- 
tions - by the emission from rather small, dense loops. 
Stellar coronae always appear as point sources. The only 
way to infer structural information in these unresolved 
point sources has been via eclipse studies in suitably 
chosen systems where one tries to constrain the emitting 
plasma volume from the observed light curve. Another 
method to infer structure in spatially unresolved data are 
spectroscopic measurements of density. The emissivity 
of plasma in coronal equilibrium in carefully selected 
lines does depend on density. Some lines may be present 
in low-density plasmas and disappear in high-density 
plasmas such as the forbidden lines in He-like triplets, 
while other lines may appear in high-density plasmas 
and be absent in low-density plasmas (such as lines 
formed following excitations from excited levels). With 
the high-resolution spectroscopic facilities onboard 
Chandra it is possible to carry out such studies for a 
wide range of X-ray sources. The purpose of this paper 
is to present and discuss some key density diagnostics 
available in the high-resolution grating spectra obtained 
with Chandra. We will specifically discuss the spectra 
obtained with the Low Energy Transmission Grating 
Spectrometer (LETGS) for the stars Capella and Procyon. 



Both Capella and Procyon are known to be rel- 
ative steady and strong X-ray sources; no signatures 
of flares from these stars have ever been reported in 
the literature. Both Capella and Procyon are rather 
close to the Sun at distances of 13 pc and 3.5 pc 
(Tab. [I]), so that effects of interstellar absorption are 
very small. Both of them have been observed with vir- 
tually all X-ray satellites flown so far. Capella was first 
detected as an X-ra y source by |Catura et al. (1975) 



and confirmed by [Mewe et al. (1975) , Procyon by 



Schmitt et al. (1985). The best coronal spectra of Capella 
were obtained with the Einstein Observatory FPCS 



and OGS QVedder et al. (1983)[ |Mewe et al. (1982)1) 



the EXOSAT trans mission grating (|Mewe et al. (1986) 
1989)[) and EUVE flDupree et al. (1993) 



Lcmcn et al. 



Schrijver et al. (1995) ), while high-spectral resolution 



Table 1. Properties of Procyon and Capella: mass M, 
radius R, effective temperature T e g, log g and the limb 
darkening coefficient e 



Procyon 



Capella 



d/pc 
M/Mq 
Rj Rq 
T cff /K 

log g 
Spectr. type 



3.5 

1.7±0.1 4 
2.06 ± 0.03 4 
6530 ± 90 2 
4.05 ± 0.04 2 
F5 IV-V 

0.724 1 



13 

2.56 ± 0.04 3 
9.2±0.4 3 
5700 ± 100 3 
2.6±0.2 5 
Ab: Gl III 
(Aa: G8/K0 III) 
0.83 1 



References: 



Di'az-Cordovcs ct al. (1995)| 



Fuhrman et al. (1997 



Hummel ct al. ( 1994)1 , 
Irwin et al. ( 1992)|, 



Kelch et al. (1978) 



Mewe et al. (1986)| and |Lemen et al. (1989)| using EX 
OS AT transmission gratin g data and |Drake et al. 1995 



and |Schrijver et al. (1995)| using EUVE data. No te that 
[Schmitt ct al. (1996b) and [Schrijver et al. (1995) investi- 
gated the coronal density of Procyon using a variety of 
density sensitive lines from Fex to Fexiv in the EUV 
range and found Procyon's coronal density consistent with 
that of solar active region densities. 

The plan of our paper is as follows: We first briefly re- 
view the atomic physics of He-like ions as applicable to 
solar (and our stellar) X-ray spectra. We briefly describe 
the spectrometer used to obtain our data, and discuss in 
quite some detail the specific procedures used in the data 
analysis, since we plan to use these methods in all our 
subsequent work on Chandra and XMM-Newton spectra. 
We then proceed to analyze the extracted spectra and de- 
scribe in detail how we dealt with the special problem of 
line blending with higher dispersion orders. Before pre- 
senting our results we estimate the formation tempera- 
tures of the lines, the influence of the stellar radiation 
field and the influence of optical depth effects followed by 
detailed interpretation. The results will then be compared 
to measurements of the Sun and we close with our conclu- 



2. Atomic physics of He-like ions 



The 
like 



theory 
triplets 



of the 
has 



atomic physics of helium- 
been extensively described 



spectral data for Procyon have been presented by 



in the literature ( 


Gabriel & Jordan (1969) 


Blumcnthal et al. (1972), 


Mewe & Schrijver (1978) 


Pradhan ct al. (1981)1 


Pradhan & Shull (1981) 


Pradhan (1982), |Pradhan (1985) 


, and recently 
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Porquet fc Dubau (2000] and [Mewe et al. (2000a)| ). 
Basically, the excited states (ls2Z) split up into the 
terms 2 1 P, 2 3 P,2 1 S, and 2 3 S, out of which the levels 
with J ^ decay to the ground state l x 5 through the 
resonance line (abbreviated by r), the intercombination 
line (i) and the forbidden line (f), respectively; the latter 
two lines involve spin changes and therefore violate the 
selection rules for electric dipole radiation. Although the 
radiative transition rate for the forbidden line is quite 
small, in a low-density plasma collisional depopulation 
processes are so rare, that the excited 2 3 S state does 
decay radiatively. In a high-density plasma collisional 
deexcitations dominate and hence the forbidden line 
disappears. Complications arise from other competing 
processes populating and depopulating the 3 P and 
3 S levels. These are in particular radiative transitions 
induced by the underlying photospheric stellar radiation 
field (discussed in Sect. |5.2| ) as well as ionization and 
recombination processes from the Li-like and H-like ions. 

It is customary to describe the measured line intensi- 
ties r, i and f in terms of the ratios 



Rot 



and 



Gobs 



f 



(1) 



In this paper we use the notation used in e., 



Pradhan fc Shull (1981) . In other contexts (e.g. calculat- 
ing ions of higher ionization stages) another notation la- 
beling r as w, f as z and i as x+y is often used pointing 
out that the intercombination line separates via M2 and 
El transitions into two components. Theory describes the 
density sensitivity of the R ratio by the functional depen- 
dence 



R(N e ) 



Bo 



1 



1 + 0/0 c + N c /N c 



(2) 



Densities are inferred from equating R(N e ) with P b s - For 
convenience, we follow h ere the expressions as given by 
Blumcnthal et al. (1972) . The low-density limit Rq, which 
applies for A c = and = 0, is given by 



Rq 



1 + F 
B 



1. 



(3) 



with the radiative branching parameter B averaged over 
the two components of the intercombination line. The pa- 
rameter N c , the so-called critical density (cf., Tab. |J), 
above which the observed line ratio is density-sensitive, 
is given by 



Nr. 



A(2 3 S 1 ^ 1%) 
(l + P)G(2 3 Si -> 2 3 P) 



(4) 



where A(2 3 Si — > l 1 ^) is the radiative transition proba- 
bility of the forbidden line, G(2 3 Si -> 2 3 P) is the elec- 
tron collisional rate coefficient for the transition from 2 3 S' 
to 2 3 P, and F is the ratio of the collision rates from the 
ground to the levels 2 3 S" and 2 3 P, respectively (including 



Table 2. Atomic Parameters for He-like triplets. T m is 
the peak line formation temperature (MEKAL), Ro is the 
low-density limit and N c the density were R falls to half 
its low-density value. 



ion 


T m /MK 


Ro 


iV c /(10 10 cm- 3 ) 


Cv 


1 


10.6 


0.051 


Nvi 


1.4 


4.9 


0.45 


Ovn 


2.0 


3.495 


3.00 



cascade effects from upper levels through radiative tran- 
sitions after excitation or recombination) . Finally, the pa- 
rameters <fi (the radiative absorption rate from 2 3 5* to 2 3 P 
induced by an external radiation field) and 4> c describe the 
additional possible influence of the stellar radi ation field 
on the depopulation of the 3 S state (cf., Sect. 5.2). The 
latter is insignificant especially for the He-like triplets of 
higher ionization stages than Cv. Values for R and N c 
used in this paper are listed in Tab. p[ Obviously, for very 
large densities or large radiation fields the forbidden line 
and hence R disappears. The ratio 



G = (i + f)/r, 



(5) 



a measure of the relative strength of the resonance 
line, depends on t he electron temperature, T (e.g., 
Mewe et al. (2000a) ). G can be used to derive tempera- 
tures for the lines used for density diagnostics that can be 
compared to the maximum formation temperature T m . 

3. Instrument description 

The Low Energy Transmission Grating Spectrometer 
(LETGS) on board the Chandra Observatory is a diffrac- 
tion grating spectrometer covering the wavelength range 
between 2-175 A (0.07-6 keV). 540 individual grating 
elements are mounted onto a toroidal ring structure. Each 
of the elements consists of a freestanding gold grating with 
1/im grating period. The fine gold wires are held by two 
different support structures, a linear grid with 25.4 /im 
and a coarse triangular mesh with 2 mm spacing. The 
whole grating ring can be inserted into the convergent 
beam just behind the High Resolution Mirror Assembly 
(HRMA) thereby dispersing the light of any X-ray source 
in the field of view into its spectrum. The efficiency of the 
grating spectrometer is of the order of 10% on average but 
is enhanced by a factor of two around 2 keV due to partial 
transparency effects; a more detailed description of the in- 



strument is presented by Predehl et al. (1997). Both sides 
of the spectrum are recorded with a microchannel plate 
detector (HRCS), placed behind the transmission grating. 
In contrast to CCD based detectors, the HRCS detector 
provides essentially no intrinsic energy resolution, the en- 
ergy information for individually recorded events is solely 
contained in the events' spatial location. 
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0.020 



0.01 



< 0.000 

-0.01 
-0.020 



Table 3. List of the data sets with start and stop times 
and durations used in this paper. 







>0 
A 



30 

A 



40 



50 



Fig. 1. Comparison of fit results from left side (A_) and 
right side (A+) of the spectrum. 



4. Observations and data analysis 

4--1- Observations 

The data described and analyzed in this paper were 
gathered during the calibration phase of the Chandra 
LETGS. Part of the da ta has already been presented by 
Brinkman ct al. (2000) . The individual observation inter- 
vals used in our analysis are listed in Tab. 0. As is clear 



from Tab. g, the observations extend over two days for 
Procyon and over two months for Capella. Our analysis 
refers to the mean properties of Procyon and Capella; 
the subject of possible time variations is not subject of 
this paper. The total on source integration times are 
218.54 ksec for Capella and 140.75 ksec for Procyon. Since 
the sources are among the strongest coronal X-ray sources 
(cf., Hiinsch ct al. (1998a, b)), it is also evident that these 
spectra belong to the highest quality coronal X-ray spec- 
tra reasonably obtainable with the Chandra LETGS. 

4-2. Data extraction 

All the HRCS datasets analyzed in this paper (see Tab. |J) 
were processed using the standard pipeline processing. 
The incoming X-rays are diffracted by the grating, dispers- 
ing the different energy photons to different detector posi- 
tions along the dispersion direction. Therefore the spectral 
information obtained with the LETGS has to be extracted 
spatially. The pulse heights with which the HRC detector 
records registered events contains some very modest en- 
ergy information, which was, however, not used. At each 
wavelength, the photons have to be integrated in cross- 
dispersion direction. Because of the spectrograph's astig- 
matism, the width of the spectral trace in cross-dispersion 
direction is wavelength dependent. For wavelengths below 



Obs. 


time 


Observation date [UT] 


ID. 


[ksec] 


start 


end 


Capella 








62435 


22.34 


09-06-1999 00:35:40 


09-06-1999 06:48:01 


01167 


15.36 


09-09-1999 13:10:06 


09-09-1999 17:26:08 


01244 


12.37 


09-09-1999 17:42:27 


09-09-1999 21:08:36 


62410 


11.33 


09-09-1999 23:43:57 


09-10-1999 02:52:48 


01246 


15.00 


09-10-1999 03:06:06 


09-10-1999 07:16:08 


62422 


11.68 


09-12-1999 18:26:42 


09-12-1999 21:41:20 


62423 


14.80 


09-12-1999 23:37:44 


09-13-1999 03:44:28 


01420 


30.30 


10-29-1999 22:49:29 


10-30-1999 07:14:27 


01248 


85.36 


11-09-1999 13:42:24 


11-10-1999 13:25:05 


total: 


218.54 






Procyon 






00063 


70.39 


11-06-1999 21:24:31 


11-07-1999 16:57:38 


01461 


70.36 


11-07-1999 17:04:55 


11-08-1999 12:37:39 


total: 


140.75 







75 A we choose a 3.6" wide extraction box around the 
spectral trace which includes almost all of the source signal 
while keeping the background level low. For wavelengths 
greater than 75 A the extraction box widens in cross dis- 
persion direction to 8.4" at 175 A. In addition the photons 
in four background regions of identical shape to that of the 
object extraction region where selected. These regions are 
displaced 12" and 24" above and below the source extrac- 
tion box in order to check for any spatial variation of the 
background in cross-dispersion direction. 



4-3. Symmetry of the grating spectrum 

For the purpose of verifying the symmetry of the dis- 
persed spectra, we determined empirical wavelengths of 
ten strong emission lines recorded on both sides of the 
middle detector plate (i.e. A < 50 A). For the example 
of Capella we plot the wavelength difference (A+ — A-) 
against the mean wavelength ^(A+ + A_) in Fig. [j]. As 
can be seen from Fig. |l|, the wavelength values obtained 
from the right and left side agree to within less than the 
instrumental resolution of ~0.06 A (cf. ct- values in Tab. 
Both spectra can thus be added for the analysis in order 
to increase the SNR of the data. In the following we will 
always consider co-added spectra. 



.4- Method of analysis 
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4.4.1. Treatment of background 

Before determining the flux of any emission line, the spec- 
tral background must be modeled. This background con- 
sists of essentially two components: instrumental back- 
ground and source background, which may consist of con- 
tinuum radiation and weak lines. The instrumental back- 
ground is also present on those parts of the microchannel 
plate not illuminated by X-rays from the target source. We 
checked for any variation of the instrumental background 
perpendicular to the dispersion direction and found none. 
We then modeled this instrumental background bg along 
the dispersion direction by a low-order polynomial and 
found acceptable fits for almost the whole spectrum. A 
far more difficult task is the determination of a reliable 
source background sbg. The source background is by def- 
inition present only on the extracted spectral trace; it is 
formed by continuum radiation from the source, and pos- 
sibly by weak, unrecognized spectral lines and/or higher 
order contamination. 

For the determination of the fluxes of individual emission 
lines, the source background is required only in the vicin- 
ity of the line(s) under consideration and is approximated 
by a single number. Our numerical procedures can, how- 
ever, scope with arbitrary source background models. In 
order to estimate this source background, we adopted the 
following procedure. We first subtracted the instrumental 
background from the spectral trace on the source, and cal- 
culated the median of the thus obtained count values for 
all bins within the investigated (small) part of the spec- 
trum. The median is a statistically robust estimate of the 
background (which we assume to be flat over the consid- 
ered part of the spectrum) as long as more than 50 % of 
all bins belong to the source background, i.e. the spectrum 
contains not too many lines. For each bin i we thus obtain 
a background value of sbg + bgi in the studied wavelength 
range. 

4.4.2. Fit procedure 

The emission line spectra are fitted with a maxi- 
mum likelihood technique similar to that used by 
Schmitt et al. (1996a) for their EUVE spectra. The spec- 



trum is assumed to consist of the background b and of M 
discrete emission lines. The treatment of the background 
b is described in Sect 



Let the observed spectrum be given on a grid of N bins 
with wavelength values Ai . . . Ajy- The number of expected 
counts Ci in the i th bin can then be calculated as 



M 



sbg + bg. t + a j g l j , 



(7) 



where a,j is the total number of counts of line j and gij is 
the value of the profile function for line j in bin i. In or- 
der to compare the modeled spectrum Cj with a measured 
spectrum rii, we assume rii to be a Poisson realization of 
Cj. The total probability of the observations n\ ... njy is 
then given by 



P{ni 



n N ) 



N 

n< 

i=l 



(8) 



and the likelihood function L(aj ,Xj,aj,j — 1 . . . N) is de- 
fined as 



JY 



C = —2 InP = — 2 Ci + rti lncj) + const . 



(9) 



The best-fit values of the parameters a,j, Xj and <jj (j = 
1 . . . N) are determined by finding extremal values of C 
through: 



dhiC 

da 3 



= 0. 



dlnC 
~dX~ 



= 0. 



d\nC 

da j 



= 0. 



(10) 



The physical meaning of the fit parameters cij, Xj and Uj 
is quite different. Xj and o~j are in principle fixed by the 
wavelengths of the considered lines and the instrumental 
resolution, complications arise from possible wavelength 
calibration errors and line blends. On the other hand, 
the amplitudes aj, proportional to the line flux, are the 
genuine interesting parameters. Analysis of the likelihood 
equations shows that for the amplitude aj a fixed point 
equation can be derived for each line j with Xj and aj 
assumed to be given: 



N 
i=l 



(11) 



Eq. 



4.4.1 



11 can be efficiently solved by iteration. In order 

and line- 



Each line j is assumed to to find °P timum values for the wavelengths Xj 



be represented by a normalized profile Qj(X; Xj, aj), e.g. a 
Gaussian profile 



«/,(A:A,.rr,j 



1 



( A — A ) 



27^7, 



(6) 



with the dispersion aj and the central wavelength Xj . The 
assumption of a Gaussian line profile is of course arbi- 
trary. Strong isolated emission lines (like Ovin Ly a ) can 
be fitted quite well with such a model, and other analyt- 
ical models can easily be incorporated into our scheme. 



widths aj , we seek minimal values of C by ordinary min- 
imization procedures. In this process the wavelengths Aj 
can vary either freely, or - if the wavelengths of the lines 
to be fitted are all known - the wavelength differences 
between the individual lines in a multiplet can be kept 
fixed in order to account for possible shifts of the overall 
wavelength scale. Similarly, the line widths aj can either 
vary freely or be fixed, and in such a way blended lines 
can be described. In this fashion an optimal value for C is 
obtained, and the parameters aj, Xj (if fitted) and aj (if 
fitted) represent our best fit measurements of these values. 
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Measurement errors are determined by assuming the like- 
lihood curve C{a,j) to be parabolic and finding the value 
of Aa,j where £(<Xj ± Aaj) = £(dj) + d£ resulting in 

Aa 7 - 



Table 5. Fit result in the 13 -14 A range for Capella. 
The assumed source background is 5792 cts/A 



^ = ^j2A£/£"- with £^ being the second derivative 
of £ with respect to aj known from Eqs. ^| and [t]. We 
choose A£ — 1 which yields formal 1 a errors. 
Similarly the errors for Aj and <jj are calculated. We de- 
termine the errors in dA,, and do~j numerically from £(Xj) 
and £{<Jj), respectively with A£ — 1, i.e. errors being given 
within 68.3%. Thus formally we treat all other parameters 
as "uninteresting". 



4-5. Extracted spectra and measured line ratios 



With the procedure described in Sect. 4.4 we analyzed the 
spectra of the He-like triplets of O VII, N vi, and C v for 
the two stars Capella and Procyon. In Figs. || and || the 
measured spectra are plotted with bold line-dotted lines 
and the best-fit model curve is indicated by a thin solid 
line. The total background, i.e. the instrumental back- 
ground and the assumed source continuum background, 
is shown with a dotted line. The derived best-fit param- 
eters are given in Tab. |4j, where we list, for both Capella 
and Procyon, the derived empirical wavelengths, the line 
widths and the line strengths (in counts). For reference 
purposes we also list the source background values used. 
Extrapolating the instrumental background onto the spec- 
tral trace shows that Procyon's source background must 
be quite small as expected for a low coronal temperature 
X-ray source. It will therefore be neglected for the purpose 
of line flux modeling. In Tab. |4| we also list the ratios be- 
tween forbidden and intercombination line, and that of 
the sum of intercombination and forbidden line to the 
resonance line, which are needed for subsequent analysis. 
Since no significant rotational or orbital line broadening 
is expected given the even high spectral resolution of our 
Chandra data, the line- width a was kept fixed for all mod- 
els assuming that this value described the instrumental 
resolution. Starting values for the wavelengths were taken 
from |Mewe et al. (1985)| . 



The O vn triplet has the best signal to noise ratio in both 
stars and is unaffected by any significant contamination 
from higher orders as can be seen from Fig. || and Fig. || 
in the top panel. As far as the N vi triplet is concerned, it 
is not as isolated as the O vn triplet and requires further 
analysis. In both Procyon (Fig. ||) and Capella (Fig. ||) 
additional lines appear. At 28.44 A a line attributed to 
C VI is evident for both Capella and Procyon. In order to 
minimize any possible cross talk, this line was included 
in the fit. The line at 30 A seen in the Capella spectrum 
is interpreted as the second order of the strong FeXVH 
line at 15.013 A; Fexvn is strong in Capella, but essen- 
tially absent in Procyon. A somewhat strange feature can 
be seen only in the Nvi spectrum of Procyon at about 



Line 


A [A] 


a [A] 


A [cts] 


Ne IX 


13.46 ± 0.003 


0.025 


4373.9 ± 82.92 


Fexix 


13.53 ±0.001 


0.025 


4683.4 ± 84.52 


Ne ix 


13.71 ±0.004 


0.033 


3019.9 ± 74.50 


Fexix 


13.8 ±0.02 


0.02 


496.6 ± 62.06 


Fe xvii 


13.83 ±0.008 


0.03 


3337.6 ± 77.94 



29.3 A. This is not a line but an instrumental effect which 
is present only on the negative side of the spectrum. 

4.5.1. Analysis of the Cv triplet 

Contamination with higher order lines makes the analysis 
of the Cv triplet for Capella particularly difficult (cf., 
Fig. ||). This is clear since the Capella spectrum contains 
strong emission in the band between 13-14 A, which 
appears in third order in the 39-42 A band under con- 
sideration for the Cv triplet. Specifically, the interfering 
lines are Neix (13.44 A x 3 = 40.32 A), Fexix (13.52 A 
x 3 = 40.56 A), Neix (13.7 A x 3 = 41.1 A), Fexix 
(13.795 A x 3 = 41.385 A) and Fexvn (13.824 A x 3 = 
41.472 A). We modeled the contamination from these five 
lines by first determining their first order contributions, 
transforming the fit results to third order. In the transfor- 
mation the first order intensities are reduced by a factor 
of 14.1. This reduction factor is obtained by comparison 
of first and third order of the isolated 15.013 A line 
(Fexvn). Since the contaminating photons are rather 
energetic, we expect essentially the same reduction 
factor in the range 13-14 A. The first order fit result is 
plotted in Fig. [I] and the resulting fit values are listed 
in Tab. ^|. Technically the higher order contamination 
was treated as an additional nonconstant contribution 
to the instrumental background. In addition to the C v 
triplet, another line, Sixn at 40.91 A, appearing strong 
in Capella but weaker in Procyon, had to be modeled. 
The source background was estimated by requiring the 
lowest count bins to be adequately modeled. The median 
function could not be applied since there are too many 
lines in the considered wavelength range. The result of 
this modeling exercise is shown in Fig. [| in the bottom 
panel, where the dotted line, representing the background, 
is not constant, but heavily influenced by higher order 
lines. We emphasise that the errors listed in Tab. ^ do 
not include errors from the fits in the 13 A band. We 
point out in particular that the forbidden C v line lies 
on top of the third order Fexvn 13.824 A line, so that 
any derivation of the forbidden line flux does require an 
appropriate modeling of the third order contamination. 
This is especially difficult because the Fexvn 13.824 A 
line is blended with the Fexix 13.795 A line in first order, 
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Table 4. Measured line ratios for Capella and Procyon 



Capella 


A [A] 


a [A] 


A [cts] 


sbg [cts/ A] 


Robs —f/i 




Ovii 
r 
i 
f 


21.62 ± 0.005 
21.82 ± 0.007 
22.12 ± 0.001 


0.027 


3071.2 ± 56.0 
544.8 ±31.4 
2135.2 ± 51.1 


2997 


3.92±0.24 


0.87±0.03 


Nvi 
r 
i 
f 


28.79 ± 0.003 
29.1 ± 0.004 
29.54 ± 0.003 


0.03 


491.2 ± 31.49 
228.2 ± 26.5 
384.5 ± 29.4 


3265 


1.68±0.23 


1.25±0.14 


Cv 
r 
i 
f 

Procyon 


40.28 ± 0.002 
40.72 ± 0.008 
41.5 ± 0.005 


0.026 
0.026 
0.029 


440.7 ± 26.9 
101.3 ± 18.24 
160.2 ± 22.37 


1200 


1.58±0.36 


0.59±0.1 


Ovii 
r 
i 
f 


21.6 ± 0.006 
21.8 ± 0.004 
22.1 ± 0.003 


0.027 


731.6 ± 28.7 
203 ± 16.8 
652.4 ± 27.3 





3.21±0.3 


1.17±0.08 


Nvi 
r 
i 
f 


28.8 ± 0.003 
29.1 ± 0.006 
29.55 ± 0.005 


0.03 


200.2 ± 16.8 
77.4 ± 12.3 
97.1 ± 13.2 





1.26±0.26 


0.87±0.14 


Cv 
r 
i 
f 


40.28 ± 0.003 
40.75 ± 0.004 
41.48 ± 0.01 


0.03 
0.028 
0.046 


203.8 ± 17.0 
123.1 ± 14.2 
63.4 ± 13.5 





0.51±0.12 


0.92±0.16 



while in third order this blend is resolved. 
The strong emission in the 13 A regime is mostly due to 
iron in excitation stages Fexvil and Fexix. Since there is 
no significant emission from iron in these high excitation 
stages in the spectrum of Procyon, there is no blending 
with third order lines in this case. From Fig. || it can be 
seen that the modeling is straightforward for Procyon. 



5. Results and interpretation 

Inspection of our Chandra spectra of Capella and Procyon 
shows the O VII and N vi forbidden lines to be present in 
both cases. The C v forbidden line is very weak but present 
for Procyon, it is also present in Capella despite being 
significantly contaminated by third order radiation. The 
observed f/i-ratios for O vn are very close to the expected 
low-density limit i?o for both Procyon and Capella, while 
the measured f/i-ratios for Nvi and Cv are definitely be- 
low the respective values of Rq in both cases. Before giving 
a quantitative interpretation of the observed line ratios we 
must first consider the formation temperature of the stud- 
ied He-like lines. Especially for the C V triplet the stellar 



radiation field can contribute significantly to the depop- 
ulation of the atomic level from which the forbidden line 
originates. Low forbidden line intensities can thus also in- 
dicate large radiation fields and not necessarily high den- 
sities. 



5.1. Formation temperature of lines 

For the purpose of density diagnostics it is customary to 
assume that all of the emission is produced at a single 
temperature T m (cf., Tab. ||), which corresponds to the 
peak of the contribution function of the considered line, 
the so-called formation temperature. It should be kept in 
mind, however, that He-like ions are present over a rel- 
atively broad temperature range, and therefore this as- 
sumption might be poor if steep emission measure gradi- 
ents are present. The measured value of G is also a tem- 
perature diagnostics, although optical depth effects in the 



resonance line may contribute (cf., Acton (1978)); the rel- 
evance of opt ical depth effects to our results is analyzed 
in section 5.4 . In Tab. || we show the temperatures T(G) 
determined from converting the observed G hs values into 
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Fig. 2. Spectrum (line-dotted) and fitted (solid) curve for 
the Ovil, Nvi, and Cv triplet for Procyon. The dotted 
line represents the total background. The binsize is 0.02 A 
for Ovil and 0.03 A for Cv and Nvi. 



Fig. 3. Spectrum (line-dotted) and fitted (solid) curve for 
the O vn, N vi, and C v triplet for Capella. The dotted line 
represents the total background. The binsize is 0.02 A for 
Ovn and 0.03 A for Cv and Nvi. 



Blending at 1 3 A 



temperatures according to Mewe et al. (2000a) using an 
electron density of n e = 5 • 10 9 cm -3 . We also determined 
temperatures by studying the line ratio between the ob- 
served Ly a lines of Ovm, Nvn, and Cvi and the r-lines 
of Ovil, Nvi, and Cv respectively by assuming isother- 
mal plasma emission. We assume plasma emissivities 



as calculated in th e codes ME KAL flMcwe et al. (1985) 
Mewe et al. (1995)[ ) and SPEX flKaastra et al. 1996| ). The 



thus obtained fluxes are multiplied with effective areas 
(cf., Tab. 0) before comparison with the measured ratios. 



as T(H-He) 
% the temperatures T(G) 



The results are listed in Tab. 

As can be seen from Tab. 
and T(H-He) do not agree. This is not surprising given 
the fact that we are likely dealing with a temperature dis- 
tribution. For Procyon the temperatures T(H-He) agree 
quite well with T m , while T(G) agrees with T m only for 
Nvi. For Capella the respective temperatures T(H-Hc) 
always exceed those found for Procyon, while the T(G) 
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Fig. 4. First order spectrum responsible for third order 
contamination of the Cv triplet of Capella; shown are 
the data (dash-dotted histogram) as well as the best fit 
(solid line). Note that the feature near 13.8 A is actually a 
line blend (consisting of Fexix and Fexvn contributions) 
resolved in third order. 
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Table 6. Temperature diagnostics from ratio of Ly a /r 
and from the ratio G=(i±f)/r. Values for r and G are 
taken from Tab. || (the G values listed are corrected for 
effective areas; cf., Tab. 0). For comparison the peak line 
formation temperature T m is also given. T(G) is calculated 



according to Mewe et al. (2000a) 



Ion 


A 


Procyon 


Capella 




[A] 


[cts] 


[cts] 


Cvi (cts) 


33.75 


697.5 ± 28.1 


2151.1 ± 51.3 


C v (cts) 


40.27 


203.8 ± 17.0 


440.7 ± 26.9 


Cvi/Cv 




3.42 ± 0.32 


4.88 ± 0.41 


T(H-He)/MK 




1.14 ± 0.03 


1.27 ± 0.03 


Gobs 




1.27 ± 0.21 


0.81 ± 0.14 


T(G) /MK 




0.32 ± 0.15 


0.98 ± 0.43 


T m /MK 




1.0 


1.0 


Nvn (cts) 


24.8 


206.9 ± 16.95 


2280.3 ± 53.6 


Nvi (cts) 


28.79 


200.2 ± 16.8 


491.2 ± 31.49 


Nvii/Nvi 




1.03 ± 0.12 


4.64 ± 0.41 


T(H-He)/MK 




1.5 ±0.05 


2.5 ±0.08 


Gobs 




0.93 ± 0.16 


1.33 ± 0.15 


T(G)/MK 




1.25 ± 0.60 


0.47 ± 0.17 


Tm/MK 




1.4 


1.4 


Oviii (cts) 


18.97 


673.2 ± 27.6 


14676.8 ± 124.3 


Ovu (cts) 


21.6 


731.6 ± 28.7 


3071.2 ± 56 


viii/O vn 




0.92 ± 0.05 


4.78 ± 0.13 


T(H-He)/MK 




2.13 ±0.03 


3.4 ±0.03 


Gobs 




1.21 ± 0.08 


0.90 ± 0.03 


T(G) /MK 




1.0 ± 0.16 


2.01 ± 0.16 


Tm/MK 




2.0 


2.0 



Table 7. Effective areas of the detector at the designated 
wavelengths. The values are taken from In-Flight Calibra- 
tion by Dcron Pease (9 March 2000). 



element 


A / A 




A off / cm 2 








A+ 


A_ 


total 


Cvi 


33.76 


5.09 


5.09 


10.18 


Cv (r) 


40.28 


2.34 


2.34 


4.68 


Cv (i) 


40.74 


1.65 


1.65 


3.30 


Cv (f) 


41.5 


1.76 


1.76 


3.52 


Nvn 


24.8 


7.76 


7.75 


15.51 


Nvi (r) 


28.8 


6.34 


6.33 


12.67 


Nvi (i) 


29.1 


6.16 


6.15 


12.31 


Nvi (f) 


29.54 


5.80 


5.79 


11.59 


O vin 


18.98 


10.08 


10.70 


20.78 


Ovu (r) 


21.62 


6.47 


7.05 


13.52 


O vn (i) 


21.82 


6.36 


6.94 


13.30 


O vn (f ) 


22.12 


6.22 


6.81 


13.03 



temperature derived from N vi is below that found for 
Procyon. We tentatively conclude that the observed Nvi 
and O vn emission in Capella has significant contributions 
from plasma at temperatures away from the peak in the 
line emissivity curve, while for Procyon the emission ap- 
pears to come from rather close to the line emissivity peak. 

5.2. Influence of the stellar radiation field 

The observed C v line ratios are particularly interest- 
ing. While the observed value for Capella agrees well 
with the solar observations (cf., Austin et al. (1966) 



Freeman & Jones (1970)), the observed value for Procyon 
is rather small. The important point to keep in mind in 
this context is that Procyon and Capella are stars with 
different properties compared to the Sun. Procyon is of 
spectral type F5V-IV (cf., Tab. |l|) with an effective tem- 
perature of 6500 K, Capella is a spectroscopic binary, the 
components of which are of spectral type Gl and G8; 
occasionally the brighter component is also classified as 
F9. Strictly speaking, what really matters is the effec- 
tive temperature of the radiation field at the wavelength 
corresponding to the energy difference between forbidden 
and intercombination line levels, i.e., 2272 A for Cv (cf., 
Tab. |) for the transition 2 3 P -> 2 3 S. We investigated 
the stellar surface radiation fluxes from measurements ob- 
tained with the International Ultraviolet Explorer satellite 
(IUE) . We first determined continuum fluxes from archival 
IUE data, and converted these fluxes into intensities using 
the expression 



h=F, 



d 2 



1 



B? 2tt(1-§)' 



(12) 



where e is the (linear) limb darkening coef ficient. It is de- 
termined from Diaz-Cordoves et al. (1995) using log g and 
T cf f listed in Tab. |l|; d and R are distance and radius of 
the sample stars as listed in Tab. |l|. For the case of the 
binary system Capella we assumed the worst case and at- 
tributed all emission to the Ab component, i.e., the star 
with the higher radiation temperature. The other compo- 
nent Aa (G8 /K0 III) is about equally bright in the corona 
as Gl III (cf, |Linsky et al. (1998) ), but much weaker in 
the chromosphere and transition region. From the thus 
obtained value for the intensity a radiation temperature 
T ra d can be inferred from the appropriate Planck curve, 
expressed as 



W 



8tt/i^ 3 



exp 



hv 



(13) 



where W is the dilution factor of the radiation field (we 
take W = \ close to the stellar surface). The results 
derived from these data are listed in Tab. 0. 
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Table 8. Investigation of the influence of the stellar radi- 
ation field. Measured fluxes from the IUE satellite F\ are 
converted to intensity I\ taking into account limb dark- 
ening effects using e from Tab. |. 



Table 9. Results of density diagnostics 







Cv 


Nvi 


O vii 




Af-i/A 


2272 


1900 


1630 




Ipot/eV 


392.1 


552.1 


739.3 






Capella 






iV(io 


— 11 ergs ^ 
cm-^ sA' 


2 ±0.5 


1.2 ±0.5 


0.25 ±0.05 


/a/(10 4 


ergs \ 
cm-^ s A strd ' 


1.98 ±0.5 


1.19 ±0.3 


0.25 ±0.5 




r rad /K 


4585 ± 100 


4976 ± 150 


5029 ± 50 




<t>/<t>c 


2.54 ±0.86 


0.2 ±0.1 


0.003 










±0.0005 






Procyon 




iV(io 


— 11 ergs \ 
cm^ sA' 


15 ±4 


7±3 


0.6 ±0.3 


/a/(10 4 


ergs \ 
cm^ s A strd ' 


21 ±5 


9.8 ±4 


0.84 ±0.4 




T r ad/K 


5532 ± 150 


5778 ± 200 


5406 ± 300 




4>/<f>c 


26.67 ±9.3 


1.58 ±0.84 


0.01 ±0.015 




^/(s- 1 ) 


34.6 


148 


717 



log n e /cm 3 


Capella 


Procyon 


Cv 


9.42 ±0.21 


<8.92 


Nvi 


9.86 ±0.12 


9.96 ±0.23 


Ovii 


<9.38 
(2a upper limit) 


q oq+0.4 



In order to compute the dependence of 
quired in Eq. f2, on T ra d we use the calculations by 



Bhuncnthal et al. (1972), who derive the expression 

4> _ 3(l±F)c 3 A(2 3 P -> 2 3 S) 
f c ~ Snhv 3 A{2 3 S 1 ^ l^o)"" 



(14) 



with <f> c = A(2 3 S X -> 
ities A(i — ► j), and u 



l 1 Sq)/(1 + F), transition probabil- 
/ the spectral energy density at the 



appropriate 2 3 Pi — > 2 3 S\ frequency v. The factor 3 is the 
ratio of the statistical weights of the levels 2 3 P and 2 3 S. 
F is approximated by Bhuncnthal et al. (1972) through 



F(0 



where 



H(0 



3gff(g)exp(C/4) + L2j + 2H(Q + 2 
^exp^/A) + 0.6£H(Z) + H($ + 1 ' 



C^S-y 2 3 S) 

C^S^ 2 3 py 



(15) 



(16) 



with £ = //fcT ra d and the collisional excitation 
rate coefficients C(i — > j). The ionization poten- 
tials /pot for the ions C v, N VI, and O vn are taken 
from 
tion 



Pradhan 



rate coefficients C(i 



Shull (1981)| . The collisional excita- 
j) were taken from 



Pradhan et al. (1981) for Cv and Ovii while the rate 
coefficients for Nvi were obtained by interpolation from 
these values. 

From Eq. 14 we determine the values used in Eq. for 
the calculation of the theoretical curves in Fig. |5|. They 
are also listed in Tab. || For comparison, we also give 
in Tab. H the values derived from recent calculations 



previously used by Porquet & Dubau (2000), and taking 
into account the effect of te mperature on G and R(n e ) 
( [Mewe et al. (2000a)| , see also [Porquet k Dubau (2000)[ ). 



5.3. Analysis and interpretation 

In Fig. H we show for Procyon and Capella the expected 
line ratio f/i as a function of the electron density n e in 
comparison with the observed line ratio i? bs (corrected 
for detector efficiencies) and its la error. The expected 
curves are plotted for the radiation temperature range 
estimated for the two stars as listed in Tab. ^; for Ovii 
the stellar radiation field does not significantly influence 
the f/i ratio as expected. We used the formation tem- 
peratures T(G) calculated from the G ratios as listed in 
Tab. ^, thus assuming all the emission being produced at 
a single temperature. In Tab. ^ we summarize the derived 
densities and their la errors for Procyon and Capella not 
accounting for errors in T(G). No density values could 
be determined for the C V triplet in Procyon and for the 
O vn triplet in Capella. Instead we give upper limits for 
the two cases (Tab. |^); for Capella only a 2a upper limit 
could be determined. 

Obviously, for both stars the measured line ratios are 
within (Capella) and very close to (Procyon) the low- 
density limit Rq. At any rate, the measured f/i ratio is 
larger for Capella than for Procyon, so one arrives at the 
somewhat unexpected conclusion that the coronal density 
in the active star Capella should be smaller than in the 
inactive star Procyon. 



5.4- Optical depth effects 

In all of the above analysis we assumed all triplet lines 
to be optically thin. In the following section we show this 
assumption to be consistent with our results. Let us there- 
fore assume that the optical depth in the resonance lines 
is not small. This leads to a reduction of the measured res- 
onance line flux due to radiative scattering. At line center 
the optical depth is given by the equation 



r = 1.2 10" 




(17) 



( gehrijver et al. 1994 ) with the fractional ionization is de- 



by Mewe et al. (2000a) considering a multi-level model noted by nj/n e ;, the elemental abundances by A z , the ra- 
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Fig. 5. Theoretical curves in the radiation temperature 
range given in Tab. || in comparison with the measured 
value of i?obs for the Ovn, Nvi and the Cv triplet for 
Capella and Procyon. The measured R values from Tab. [| 
are corrected for detector efficiencies. 

tio of hydrogen to electron density is nu/n e — 0.85, the 
oscillator strength / = 0.7 for all ions, the atomic number 
is M, the wavelength is measured in A, the temperature 
in K, the electron density n e in cm -3 and the mean free 
path is denoted by I. In Eq. [l7] we adopt a value of unity 
for the fractional ionization and use the solar abundances. 
We further assume T at the peak line formation and note, 
that t is rather insensitive to the precise value of T. We 
can determine - for each resonance line - that value of 
n e £ which yields an optical depth of unity. Adopting a 
maximum value of n e of 10 cm" " 3 (cf. Tab. |), we de- 
termine lower values of I of 1.310 12 cm, 5.610 12 cm and 
6.010 11 cm, for Ovn, Nvi and Cv, respectively. Assum- 
ing a geometry most suitable for resonance scattering, we 
can compute the respective emission measures of n\i 3 re- 
spectively. Comparing these emission measures with those 
derived from the measured line fluxes f\ 



EM 



Px(T) 



(18) 



with the line cooling function p\{T) and the distance d (cf. 
Tab.0), shows the former to be larger by a few orders of 
magnitude for both stars. This inconsistency shows that 
the assumption of a non-negligible optical depth is invalid 
and we conclude that optical depth effects are irrelevant 
for the analysis of He-like triplets in Procyon and Capella. 

6. Conclusions 

6.1. Comparison with the Sun 

Both Capella and Procyon are considered to be solar-like 
stars. In the solar context the He-like triplets of oxygen, 
nitrogen and carbon have been known for a long time and 
in fact the He-like ion density diagnostics have been first 
developed to interpret these solar data. Many solar obser- 
vations are available for the Ovn triplet, while only very 
few observations have been made for the Nvi and Cv 
triplets. The first observations of the Cv triplet are re- 
ported by Austin et al. (1966)], who o btained an f/i-ratio 
of 1.9, while Freeman & Jones (1970) found an f/i-ratio of 
1.0. The latter authors also observed the Nvi triplet, but 
failed to detect the intercombination line and hence de- 
duced f/i > 1.9. |Brown et al. (1986)| observed the C v and 
Nvi triplets during a flare and found a mean value of f/i 
2 with large scatter between 0.1-4 around this mean 
for Nvi and a value of 0.21 (with a scatter between 0.17- 
0.25) for Cv. In consequence, little can be said about the 
solar N vi line ratios because of the large measurement er- 
rors. The measured ratios for both Capella and Procyon 
are certainly consistent with the range of Nvi f/i values 



quoted by Brown et al. (1986), As far as C v is concerned, 



the flare data yield very low f/i ratios indicative of high 
densities. The Cv f/i ratio observed for Capella is higher, 
and that for Procyon - s till higher than that observed for 
the Brown et al. (1986) flare - we argue is due to large 



photospheric radiation fluxes. We therefore conclude that 
the layers contributing the Cv emission in Capella and 
Procyon are at lower density than those encountered in 
solar flares, and the same applies to the layers emitting 
Nvi. 

Observations of the O VII triplet in the solar 
corona are reported by Freeman & Jones (1970) 



McKenzie et al. (1978), McKenzic et al. (1982) 



of 
of 



which 
the 



refer to the 
observations 



and 


Brown et al. 


(1986)| 


Parkinson (1975) 


refer 1 



quiescent 
from 
refer to 
to active 



corona, 



all 

while some 



McKenzie et al. (1982) 



wo 



flares and those of 
regions. In addition 
points collected from various 



considered some data 
sources cited by Doyle (1980). In Fig. || we plot G vs 
R of these solar observations and compare these solar 
data with our Chandra measurements for Procyon and 
Capella. For clarity we omitted individual error bars, 
which are typically 0.1-0.3 for the solar measurements. 
As can be seen from Fig. ^|, most of the solar measure- 
ments yield R value between 3 and 4. A few values are 
significantly lower around R sa 2, with all them referring 
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Fig. 6. 

Ovn: 



Measurements of solar values for He-like ratios of 



Freeman fe Jones (1970) (dotted open box), 

quiescent (diamond), 
McKenzie et al. (1978) (dotted diamond), 
Doyle (1980)| (fille d triangle upside down) 



flare (asterix), 



Parkinson (1975): active region (plus symbol), 



Brown et al. (1986): flare (cross). 

Encircled are our measurements of Procyon (open triangle) and 
Capella (open box). 



to flares or active regions. Our Chandra measurements 
(3.28 for Procyon, 4.0 for Capella) thus fall into the 
bulk part of the solar data. Unfortunately, the solar 
data are by no means unambiguous. It is not always 
clear which regions the data refer to (except for the flare 
observations by Brown et al. (1986)] ) , and further, most 
of the observations are relatively old and have not been 
taken with high spatial resolution; they are comparable 
in some sense to our full disk observations of stellar X-ray 
sources. In addition, the instrumentation used in different 
experiments is quite different and each experiment is 
affected by its own specific difficulties. Nevertheless it 
seems fair to state that most O VII measurements referring 
to quiescent conditions are close to the low density limit, 
thus indicating that the physical properties of the O VII 
emitting layers in Capella and Procyon should not be 
that different from those in the Sun; a little puzzling in 
this context is the large scatter of the solar G data. In 
contrast, the solar flare data (R ~ 2) is much lower than 
any of the other measurements, consistent with a high 
density plasma. 



6.2. Comparison between Procyon and Capella 

Inspection of Tab. ^| suggests that np rocyon > ncapeiia 
for the densities derived from Cv, Nvi, and Ovn; there 
is definitely no evidence for np rocyon < ncapeiia- These 
findings are in contradiction to the claims made by 



Dupree et al. (1993), who argue on the basis of Fexxi line 
ratios measured with the EUVE satellite that ncapeiia ~ 
10 13 cm -3 . A recent analysis of the long wavelength 
part (A > 90 A) of the LETGS spectrum of Capella by 
Mewe et al. (2000b), however, yielded plasma densities of 
5 x 10 lz cm~ 3 at most for the plasma emitting in the 
Fe xxi lines with no evidence for any deviations from the 
low density limit. Of course, these values refer to those 
plasma layers emitting in the Fexxi line, which ought to 
be at far higher temperature (« 8 MK) than those emit- 
ting in the Cv, Nvi, and Ovn lines (1-2 MK). If we 
assume all the observed emissions to come from the same 
magnetically confined structures (loops), one expects the 
gas pressure to be more or less constant and hence the den- 
sity at lower temperatures even higher. However, plasma 
densities (for Capella) of ss 10 13 cm' 3 at temperatures of 
2 MK are clearly inconsistent with our LETGS spectra. 
Therefore the plasma emitting in the Fe XXI line is either 
at the low density limit or the Fe XXI emission comes from 
different structures than the emissions from He-like ions. 

Our temperature measurements of the line forming re- 
gions (cf., Tab. ^) are not fully conclusive, although they 
suggest that the same lines of Cv, Nvi, and Ovn respec- 
tively are formed at a somewhat higher temperature (at 
most a factor of 2) in Capella as compared to Procyon. 
As a consequence we must have pprocyon ~ Pcapeiia for 
the gas pressures in the respective line forming regions. 
Capella's peak coronal temperature is clearly much higher 
than Procyon's. If we again assume - going along with the 
solar analogy - that for both Procyon and Capella (all of) 
the X-ray emission originates from magnetically confined 
plasma loops, the loop top temperatures in Capella must 
be higher than in Procyon. The loop scaling laws then im- 
ply for the typical loop length scales L for Procyon and 

^ 11 LcapeLLa / r maT ^^i^ \3 A _ 



L Pr 



Capella 



mate yields 



-). A conservative esti- 



.Capella 



> 2, thus L 



Capella 



> 8 L 



Procyon- 



The conclusion then appears inescapable that the typical 
sizes of the magnetic structures in Procyon and Capella 
are quite different by probably at least one order of mag- 
nitude. 

Next, it is instructive to compare the total mean sur- 
face fluxes emitted in the helium- and hydrogen-like ions of 
carbon, nitrogen and oxygen for Procyon and Capella. Us- 
ing the measured count ratios and correcting for the effects 
of distance, exposure times and the surface areas of the 
stars we find for the C v, C vi, N vi, N VII, O VII, and O VIII 
lines values of (Capella)/ (Procyon) of 0.96, 1.37, 1.09, 4.9, 
1.87, and 9.71 respectively. Therefore up to those tem- 
peratures where the Cv, Cvi, and Nvi lines are formed, 
the mean surface fluxes in the two stars hardly differ at 
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all; pressures and temperatures are also quite similar. We 
therefore conclude that the physical characteristics and 
global properties (such as filling factor) of Procyon's and 
Capella's corona at a level of T ~ 10 6 K are very similar. 

It is interesting to again perform a comparison to 
the Sun. Freeman & Jones (1970) quote - for their SL801 
rocket flight on Nov 20 1969 - resonance line fluxes of 4.2 
10" 3 , 2.1 10~ 3 , and 3.1 10~ 3 erg/(cm 2 sec) for O VII, Nvi, 
and C v respectively. Correcting for the distance between 
Earth and Sun we find average surface fluxes of 200, 100, 
and 150 erg/cm 2 /sec for Ovu, Nvi, and C v respectively 
on the solar surface. Carrying out the same calculation 
for our Chandra data, we find - using the effective areas 
given m Tab. @ - mean surface fluxes in O vn of 3700 and 
1980 erg/ (cm 2 sec), in Nvi of 480 and 430 erg/ (cm 2 sec), 
and in Cv of 830 and 860 erg/(cm 2 sec) for Capella and 
Procyon respectively. We therefore recover our previous 
finding that Capella's and Procyon's surface fluxes are 
quite similar, and determine in addition that both stars 
exceed typical (?) solar values by one order of magnitude. 
Since the physical conditions of the line emitting regions 
are quite similar, we conclude that the coronal filling fac- 
tors are larger. 

Evaluating now the coronal pressure for Procyon we 
find - using the densities and temperatures derived from 
Nvi - p(Nvi) = 4.4 dyn/cm 2 , a value also very typi- 
cal for solar active regions. From the loop scaling law 
( [Rosner et al. (1978)D T max = 1.4 x 10 3 (pL)( 1 / 3 ) we de- 
duce a typical length of 2.7 x 10 8 cm using T max — 
1.5 x 10 6 K. This value must be considered as a lower limit 
to the probable length scale because of the unfortunate 
sensitive dependence of L on T max through L ~ T^ ax ; 
T m ax is likely somewhat higher than T max = 1.5 x 10 6 K 
as suggested by the temperatures derived from the O VIII 
data. At any rate, however, the conclusion is, that the 
X-ray emission originates from low-lying loops and hence 
also our Chandra data support the view that Procyon's 
corona has an appearance very similar to the Sun's, just 
like the conclusion drawn by Schmitt et al. (1996b) from 
their EUVE spectra. It then follows that the sizes of the 
magnetic loops in Procyon's corona are similar to the loops 
found in the solar corona. Adopting - for argument's sake 
- a characteristic length of Lp roc — 1 x 10 9 cm, we then 
find for Capella Lc ap eiia > 8 x 10 9 cm, and possibly even a 
significant fraction of Capella's radius. This result is some- 
what puzzling. On the Sun, loop structures, i.e., magnet- 
ically closed topologies, are occasionally found at larger 
heights, but they are always of low density and they never 
contribute significantly to the overall X-ray emission. It 
will be interesting to see whether this behavior is typical 
for active stars in general, or whether it applies only to 
Capella. After all, Capella may not be the prototypical 
active star. Its corona is not as hot as that of other stars, 
it does not produce flares and its radio emission is very 
weak. 
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